Chromosome-level genome assembly of milk thistle (Silybum marianum (L.) Gaertn.).

Silybum marianum (L.) Gaertn., commonly known as milk thistle, is a medicinal plant belonging to the Asteraceae family. This plant has been recognized for its medicinal properties for over 2,000 years. However, the genome of this plant remains largely undiscovered, having no reference genome at a chromosomal level. Here, we assembled the chromosome-level genome of S. marianum, allowing for the annotation of 53,552 genes and the identification of transposable elements comprising 58% of the genome. The genome assembly from this study showed 99.1% completeness as determined by BUSCO assessment, while the previous assembly (ASM154182v1) showed 36.7%. Functional annotation of the predicted genes showed 50,329 genes (94% of total genes) with known protein functions in public databases. Comparative genome analysis among Asteraceae plants revealed a striking conservation of collinearity between S. marianum and C. cardunculus. The genomic information generated from this study will be a valuable resource for milk thistle breeding and for use by the larger research community.

High-quality chromosome-level genome assembly of Nicotiana benthamiana.

Nicotiana benthamiana is a fundamental model organism in plant research. Recent advancements in genomic sequencing have revealed significant intraspecific genetic variations. This study addresses the pressing need for a precise genome sequence specific to its geographic origin by presenting a comprehensive genome assembly of the N. benthamiana LAB strain from the Republic of Korea (NbKLAB). We compare this assembly with the widely used NbLAB360 strain, shedding light on essential genomic differences between them. The outcome is a high-quality, chromosome-level genome assembly comprising 19 chromosomes, spanning 2,762 Mb, with an N50 of 142.6 Mb. Comparative analyses revealed notable variations, including 46,215 protein-coding genes, with an impressive 99.5% BUSCO completeness score. Furthermore, the NbKLAB assembly substantially improved the QV from 33% for NbLAB360 to 49%. This refined chromosomal genome assembly for N. benthamiana, in conjunction with comparative insights, provides a valuable resource for genomics research and molecular biology. This accomplishment forms a strong foundation for in-depth exploration into the intricacies of plant genetics and genomics, improved precision, and a comparative framework.

Pre-breeding of spontaneous Robertsonian translocations for density planting architecture by transferring Agropyron cristatum chromosome 1P into wheat.

KEY MESSAGE: We developed T1AL·1PS and T1AS·1PL Robertsonian translocations by breakage-fusion mechanism based on wheat-A. cristatum 1P(1A) substitution line with smaller leaf area, shorter plant height, and other excellent agronomic traits Agropyron cristatum, a wild relative of wheat, is a valuable germplasm resource for improving wheat genetic diversity and yield. Our previous study confirmed that the A. cristatum chromosome 1P carries alien genes that reduce plant height and leaf size in wheat. Here, we developed T1AL·1PS and T1AS·1PL Robertsonian translocations (RobTs) by breakage-fusion mechanism based on wheat-A. cristatum 1P (1A) substitution line II-3-1c. Combining molecular markers and cytological analysis, we identified 16 spontaneous RobTs from 911 F2 individuals derived from the cross of Jimai22 and II-3-1c. Fluorescence in situ hybridization (FISH) was applied to detect the fusion structures of the centromeres in wheat and A. cristatum chromosomes. Resequencing results indicated that the chromosomal junction point was located at the physical position of Triticum aestivum chromosome 1A (212.5 Mb) and A. cristatum chromosome 1P (230 Mb). Genomic in situ hybridization (GISH) in pollen mother cells showed that the produced translocation lines could form stable ring bivalent. Introducing chromosome 1PS translocation fragment into wheat significantly increased the number of fertile tillers, grain number per spike, and grain weight and reduced the flag leaf area. However, introducing chromosome 1PL translocation fragment into wheat significantly reduced flag leaf area and plant height with a negative effect on yield components. The pre-breeding of two spontaneous RobTs T1AL·1PS and T1AS·1PL was important for wheat architecture improvement.

QTL Mapping Using RIL Population.

Genetic maps are an excellent tool for the analysis of important traits, the development of which is the result of the combined expression of several genes, enabling the genomic localization of the factors determining them. Such features, characterized by a normal distribution of values, are referred to as quantitative or polygenic. The analysis of their genetic background using a chromosome map is called the mapping of quantitative traits loci (QTL). QTL analysis is a statistical method of determining the genetic association of phenotypic data (trait measurements) with genotypic data (DNA markers assigned to linkage groups).There are numerous tools developed for QTL mapping. This chapter introduces Windows QTL Cartographer with Composite Interval Mapping (CIM) method, which estimates the QTL position by combining interval mapping with multiple regression. The genotypic and phenotypic data used in the exemplary QTL mapping procedure were obtained for the recombinant inbred line (RIL) population of rye. Plant height, assessed in three seasons, was the exemplary trait under study.

Genetic Map Construction Using F2 and RIL/DH Mapping Population.

Genetic mapping is the determination of the position and relative genetic distance between genes or molecular markers in the chromosomes of a particular species. The construction of genetic maps uses data from the genotyping of the mapping population. Among the different mapping populations used, two are relatively common: the F2 and recombinant inbred lines (RILs) obtained as a result of the controlled crossing of genetically diverse parental forms (e.g., inbred lines). Also, the dihaploid (DH) population is often used in plants, but obtaining DHs in different crops, including rye, is very difficult or even impossible. Any molecular marker system can be used for genotyping. Polymorphic markers are used for linkage analysis, differentiating parental forms with segregation in the mapping population, consistent with the appropriate single-gene model. A genetic map is a great source of information on a species and can be an exquisite tool for analyzing important quantitative traits (QT).This chapter presents the procedure of genetic map construction with two different algorithms using the JoinMap5.0 program. First, the Materials section briefly informs about the mapping program, showing how to obtain a mapping population and prepare data for mapping. Finally, the Methods section describes the protocol for the mapping procedure itself.

Chromosome-scale genome assembly of bread wheat's wild relative Triticum timopheevii.

Wheat (Triticum aestivum) is one of the most important food crops with an urgent need for increase in its production to feed the growing world. Triticum timopheevii (2n = 4x = 28) is an allotetraploid wheat wild relative species containing the At and G genomes that has been exploited in many pre-breeding programmes for wheat improvement. In this study, we report the generation of a chromosome-scale reference genome assembly of T. timopheevii accession PI 94760 based on PacBio HiFi reads and chromosome conformation capture (Hi-C). The assembly comprised a total size of 9.35 Gb, featuring a contig N50 of 42.4 Mb and included the mitochondrial and plastid genome sequences. Genome annotation predicted 166,325 gene models including 70,365 genes with high confidence. DNA methylation analysis showed that the G genome had on average more methylated bases than the At genome. In summary, the T. timopheevii genome assembly provides a valuable resource for genome-informed discovery of agronomically important genes for food security.

A wild rice CSSL population facilitated identification of salt tolerance genes and rice germplasm innovation.

Salt stress is one of the major factors that limits rice production. Therefore, identification of salt-tolerant alleles from wild rice is important for rice breeding. In this study, we constructed a set of chromosome segment substitution lines (CSSLs) using wild rice as the donor parent and cultivated rice Nipponbare (Nip) as the recurrent parent. Salt tolerance germinability (STG) was evaluated, and its association with genotypes was determined using this CSSL population. We identified 17 QTLs related to STG. By integrating the transcriptome and genome data, four candidate genes were identified, including the previously reported AGO2 and WRKY53. Compared with Nip, wild rice AGO2 has a structure variation in its promoter region and the expression levels were upregulated under salt treatments; wild rice WRKY53 also has natural variation in its promoter region, and the expression levels were downregulated under salt treatments. Wild rice AGO2 and WRKY53 alleles have combined effects for improving salt tolerance at the germination stage. One CSSL line, CSSL118 that harbors these two alleles was selected. Compared with the background parent Nip, CSSL118 showed comprehensive salt tolerance and higher yield, with improved transcript levels of reactive oxygen species scavenging genes. Our results provided promising genes and germplasm resources for future rice salt tolerance breeding.

Mining genomic regions associated with agronomic and biochemical traits in quinoa through GWAS.

Quinoa (Chenopodium quinoa Willd.), an Andean crop, is a facultative halophyte food crop recognized globally for its high nutritional value and plasticity to adapt to harsh conditions. We conducted a genome-wide association study on a diverse set of quinoa germplasm accessions. These accessions were evaluated for the following agronomic and biochemical traits: days to 50% flowering (DTF), plant height (PH), panicle length (PL), stem diameter (SD), seed yield (SY), grain diameter (GD), and thousand-grain weight (TGW). These accessions underwent genotyping-by-sequencing using the DNBSeq-G400R platform. Among all evaluated traits, TGW represented maximum broad-sense heritability. Our study revealed average SNP density of ≈ 3.11 SNPs/10 kb for the whole genome, with the lowest and highest on chromosomes Cq1B and Cq9A, respectively. Principal component analysis clustered the quinoa population in three main clusters, one clearly representing lowland Chilean accessions, whereas the other two groups corresponded to germplasm from the highlands of Peru and Bolivia. In our germplasm set, we estimated linkage disequilibrium decay to be ≈ 118.5 kb. Marker-trait analyses revealed major and consistent effect associations for DTF on chromosomes 3A, 4B, 5B, 6A, 7A, 7B and 8B, with phenotypic variance explained (PVE) as high as 19.15%. Nine associations across eight chromosomes were also found for saponin content with 20% PVE by qSPN5A.1. More QTLs were identified for PL and TGW on multiple chromosomal locations. We identified putative candidate genes in the genomic regions associated with DTF and saponin content. The consistent and major-effect genomic associations can be used in fast-tracking quinoa breeding for wider adaptation across marginal environments.

[Genetic dissection of rice resistance to bacterial blight].

Bacterial blight, a major disease in rice, poses a serious impact on rice production. In this study, a doubled haploid (DH) population derived from a cross between the introduced japonica cultivar 'Maybelle' and the indica landrace 'Baiyeqiu' was used to investigate the pathogenicity of four pathogen races causing bacterial blight. The results showed that the pathogenicity of all the pathogen races exhibited continuous, transgressive distribution in the DH population. Moreover, strong correlations existed between every two pathogen races, with the correlation coefficients ranging from 0.3 to 0.6. A total of 12 quantitative trait loci (QTLs) distributed on chromosomes 1, 2, 3, 5, 6, 7, 9, and 12 were detected for rice bacterial blight, explaining 4.95% to 16.05% of the phenotype. Among these QTLs, a major QTL located in the interval RM6024-RM163 on chromosome 5 was detected in three pathogen races. In addition, the pyramiding of the positive alleles can apparently improve the rice resistance to bacterial blight. This study is of great significance for broadening the genetic resources with resistance to bacterial blight in China.

Repeated co-option of HMG-box genes for sex determination in brown algae and animals.

In many eukaryotes, genetic sex determination is not governed by XX/XY or ZW/ZZ systems but by a specialized region on the poorly studied U (female) or V (male) sex chromosomes. Previous studies have hinted at the existence of a dominant male-sex factor on the V chromosome in brown algae, a group of multicellular eukaryotes distantly related to animals and plants. The nature of this factor has remained elusive. Here, we demonstrate that an HMG-box gene acts as the male-determining factor in brown algae, mirroring the role HMG-box genes play in sex determination in animals. Over a billion-year evolutionary timeline, these lineages have independently co-opted the HMG box for male determination, representing a paradigm for evolution's ability to recurrently use the same genetic "toolkit" to accomplish similar tasks.

The complex polyploid genome architecture of sugarcane.

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

Chromosomal mapping of a major genetic locus from Agropyron cristatum chromosome 6P that influences grain number and spikelet number in wheat.

KEY MESSAGE: A novel locus on Agropyron cristatum chromosome 6P that increases grain number and spikelet number was identified in wheat-A. cristatum derivatives and across 3 years. Agropyron cristatum (2n = 4x = 28, PPPP), which has the characteristics of high yield with multiple flowers and spikelets, is a promising gene donor for wheat high-yield improvement. Identifying the genetic loci and genes that regulate yield could elucidate the genetic variations in yield-related traits and provide novel gene sources and insights for high-yield wheat breeding. In this study, cytological analysis and molecular marker analysis revealed that del10a and del31a were wheat-A. cristatum chromosome 6P deletion lines. Notably, del10a carried a segment of the full 6PS and 6PL bin (1-13), while del31a carried a segment of the full 6PS and 6PL bin (1-8). The agronomic characterization and genetic population analysis confirmed that the 6PL bin (9-13) brought about an increase in grain number per spike (average increase of 10.43 grains) and spikelet number per spike (average increase of 3.67) over the three growing seasons. Furthermore, through resequencing, a multiple grain number locus was mapped to the physical interval of 593.03-713.89 Mb on chromosome 6P of A. cristatum Z559. The RNA-seq analysis revealed the expression of 537 genes in the del10a young spike tissue, with the annotation indicating that 16 of these genes were associated with grain number and spikelet number. Finally, a total of ten A. cristatum-specific molecular markers were developed for this interval. In summary, this study presents novel genetic material that is useful for high-yield wheat breeding initiatives to meet the challenge of global food security through enhanced agricultural production.

QTL analysis of femaleness in monoecious spinach and fine mapping of a major QTL using an updated version of chromosome-scale pseudomolecules.

Although spinach is predominantly dioecious, monoecious plants with varying proportions of female and male flowers are also present. Recently, monoecious inbred lines with highly female and male conditions have been preferentially used as parents for F1-hybrids, rather than dioecious lines. Accordingly, identifying the loci for monoecism is an important issue for spinach breeding. We here used long-read sequencing and Hi-C technology to construct SOL_r2.0_pseudomolecule, a set of six pseudomolecules of spinach chromosomes (total length: 879.2 Mb; BUSCO complete 97.0%) that are longer and more genetically complete than our previous version of pseudomolecules (688.0 Mb; 81.5%). Three QTLs, qFem2.1, qFem3.1, and qFem6.1, responsible for monoecism were mapped to SOL_r2.0_pseudomolecule. qFem3.1 had the highest LOD score and corresponded to the M locus, which was previously identified as a determinant of monoecious expression, by genetic analysis of progeny from female and monoecious plants. The other QTLs were shown to modulate the ratio of female to male flowers in monoecious plants harboring a dominant allele of the M gene. Our findings will enable breeders to efficiently produce highly female- and male-monoecious parental lines for F1-hybrids by pyramiding the three QTLs. Through fine-mapping, we narrowed the candidate region for the M locus to a 19.5 kb interval containing three protein-coding genes and one long non-coding RNA gene. Among them, only RADIALIS-like-2a showed a higher expression in the reproductive organs, suggesting that it might play a role in reproductive organogenesis. However, there is no evidence that it is involved in the regulation of stamen and pistil initiation, which are directly related to the floral sex differentiation system in spinach. Given that auxin is involved in reproductive organ formation in many plant species, genes related to auxin transport/response, in addition to floral organ formation, were identified as candidates for regulators of floral sex-differentiation from qFem2.1 and qFem6.1.

Cytomolecular trends in Chamaecrista Moench (Caesalpinioideae, Leguminosae) diversification.

Chamaecrista is a Pantropical legume genus of the tribe Cassieae, which includes six other genera. In contrast to most of the other Cassieae genera, Chamaecrista shows significant variability in chromosome number (from 2n = 14 to 2n = 56), with small and morphologically similar chromosomes. Here, we performed a new cytomolecular analysis on chromosome number, genome size, and rDNA site distribution in a molecular phylogenetic perspective to interpret the karyotype trends of Chamaecrista and other two genera of Cassieae, seeking to understand their systematics and evolution. Our phylogenetic analysis revealed that Chamaecrista is monophyletic and can be divided into four major clades corresponding to the four sections of the genus. Chromosome numbers ranged from 2n = 14, 16 (section Chamaecrista) to 2n = 28 (sections Absus, Apoucouita, and Baseophyllum). The number of 5S and 35S rDNA sites varied between one and three pairs per karyotype, distributed on different chromosomes or in synteny, with no obvious phylogenetic significance. Our data allowed us to propose x = 7 as the basic chromosome number of Cassieae, which was changed by polyploidy generating x = 14 (sections Absus, Apoucouita, and Baseophyllum) and by ascending dysploidy to x = 8 (section Chamaecrista). The DNA content values supported this hypothesis, with the genomes of the putative tetraploids being larger than those of the putative diploids. We hypothesized that ascending dysploidy, polyploidy, and rDNA amplification/deamplification are the major events in the karyotypic diversification of Chamaecrista. The chromosomal marks characterized here may have cytotaxonomic potential in future studies.

Localization and quantitative distribution of a chromatin structural protein Topoisomerase II on plant chromosome using HVTEM and UHVTEM.

Topoisomerase II (TopoII) is an essential structural protein of the metaphase chromosome. It maintains the axial compaction of chromosomes during metaphase. It is localized at the axial region of chromosomes and accumulates at the centromeric region in metaphase chromosomes. However, little is known about TopoII localization and distribution in plant chromosomes, except for several publications. We used high voltage transmission electron microscopy (HVTEM) and ultra-high voltage transmission electron microscopy (UHVTEM) in conjunction with immunogold labeling and visualization techniques to detect TopoII and investigate its localization, alignment, and density on the barley chromosome at 1.4 nm scale. We found that HVTEM and UHVTEM combined with immunogold labeling is suitable for the detection of structural proteins, including a single molecule of TopoII. This is because the average size of the gold particles for TopoII visualization after silver enhancement is 8.9 ± 3.9 nm, which is well detected. We found that 31,005 TopoII molecules are distributed along the barley chromosomes in an unspecific pattern at the chromosome arms and accumulate specifically at the nucleolus organizer regions (NORs) and centromeric region. The TopoII density were 1.32-fold, 1.58-fold, and 1.36-fold at the terminal region, at the NORs, and the centromeric region, respectively. The findings of TopoII localization in this study support the multiple reported functions of TopoII in the barley metaphase chromosome.

Chromosome level genome assembly of endangered medicinal plant Anisodus tanguticus.

Anisodus tanguticus is a medicinal herb that belongs to the Anisodus genus of the Solanaceae family. This endangered herb is mainly distributed in Qinghai-Tibet Plateau. In this study, we combined the Illumina short-read, Nanopore long-read and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to de novo assemble the A. tanguticus genome. A high-quality chromosomal-level genome assembly was obtained with a genome size of 1.26 Gb and a contig N50 of 25.07 Mb. Of the draft genome sequences, 97.47% were anchored to 24 pseudochromosomes with a scaffold N50 of 51.28 Mb. In addition, 842.14 Mb of transposable elements occupying 66.70% of the genome assembly were identified and 44,252 protein-coding genes were predicted. The genome assembly of A. tanguticus will provide genetic repertoire to understand the adaptation strategy of Anisodus species in the plateau, which will further promote the conservation of endangered A. tanguticus resources.

Genome-wide association analysis identifies a consistent QTL for powdery mildew resistance on chromosome 3A in Nordic and Baltic spring wheat.

KEY MESSAGE: QPm.NOBAL-3A is an important QTL providing robust adult plant powdery mildew resistance in Nordic and Baltic spring wheat, aiding sustainable crop protection and breeding. Powdery mildew, caused by the biotrophic fungal pathogen Blumeria graminis f. sp. tritici, poses a significant threat to bread wheat (Triticum aestivum L.), one of the world's most crucial cereal crops. Enhancing cultivar resistance against this devastating disease requires a comprehensive understanding of the genetic basis of powdery mildew resistance. In this study, we performed a genome-wide association study (GWAS) using extensive field trial data from multiple environments across Estonia, Latvia, Lithuania, and Norway. The study involved a diverse panel of recent wheat cultivars and breeding lines sourced from the Baltic region and Norway. We identified a major quantitative trait locus (QTL) on chromosome 3A, designated as QPm.NOBAL-3A, which consistently conferred high resistance to powdery mildew across various environments and countries. Furthermore, the consistency of the QTL haplotype effect was validated using an independent Norwegian spring wheat panel. Subsequent greenhouse seedling inoculations with 15 representative powdery mildew isolates on a subset of the GWAS panel indicated that this QTL provides adult plant resistance and is likely of race non-specific nature. Moreover, we developed and validated KASP markers for QPm.NOBAL-3A tailored for use in breeding. These findings provide a critical foundation for marker-assisted selection in breeding programs aimed at pyramiding resistance QTL/genes to achieve durable and broad-spectrum resistance against powdery mildew.

Investigation of chromosomal genetic characteristics and identification of structural variation in the offspring of hexaploid triticale×hexaploid wheat.

Hexaploid triticale is an important genetic resource for genetic improvement of common wheat, which can broaden the genetic basis of wheat. In order to lay a foundation for the subsequent research and utilization of triticale germplasm materials, the chromosomal genetic characteristics of cross and backcross offspring of hexaploid triticale×hexaploid wheat were investigated in the process of transferring rye chromatin from hexaploid triticale to hexaploid wheat. Hybrid and backcross combinations were prepared with hexaploid triticale 16yin171 as the maternal parent and hexaploid wheat Chuanmai62 as the paternal parent. The chromosomes in root tip cells of F1, BC1F1 and BC1F2 plants were traced and identified non-denaturing florescence in situ hybridization (ND-FISH). The results indicated that the backcross setting rate of hybrid F1 was 2.61%. The transmission frequency of 2R chromosome was the highest in BC1F1 plants while the transmissibility of rye chromosome in BC1F2 plant was 6R>4R>2R, and the 5B-7B wheat translocation in BC1F2 plants showed severe segregation. A total of 24 structural variant chromosomes were observed both in BC1F1 and BC1F2 plants, including chromosome fragments, isochromosomes, translocations, and dicentric chromosomes. In addition, the seed length and 1000-grain weight of some BC1F2 plants were better than that of the hexaploid wheat parent Chuanmai 62. Therefore, multiple backcrosses should be adopted as far as possible to make the rapid recovery of group D chromosomes, ensuring the recovery of fertility in offspring, when hexaploid tritriale is used as a bridge to introduce rye genetic material into common wheat. At the same time, the potential application value of chromosomal structural variation materials should be also concerned.

Evolution of a plant sex chromosome driven by expanding pericentromeric recombination suppression.

Recombination suppression around sex-determining gene(s) is a key step in evolution of sex chromosomes, but it is not well understood how it evolves. Recently evolved sex-linked regions offer an opportunity to understand the mechanisms of recombination cessation. This paper analyses such a region on Silene latifolia (Caryophyllaceae) sex chromosomes, where recombination was suppressed in the last 120 thousand years ("stratum 3"). Locating the boundaries of the stratum 3 in S. latifolia genome sequence revealed that this region is far larger than assumed previously-it is about 14 Mb long and includes 202 annotated genes. A gradient of X:Y divergence detected in the stratum 3, with divergence increasing proximally, indicates gradual recombination cessation, possibly caused by expansion of pericentromeric recombination suppression (PRS) into the pseudoautosomal region. Expansion of PRS was also the likely cause for the formation of the older stratum 2 on S. latifolia sex chromosomes. The role of PRS in sex chromosome evolution has been underappreciated, but it may be a significant factor, especially in the species with large chromosomes where PRS is often extensive.

QTL epistasis plays a role of homeostasis on heading date in rice.

If there was no gene interaction, the gene aggregation effect would increase infinitely with the increase of gene number. Epistasis avoids the endless accumulation of gene effects, playing a role of homeostasis. To confirm the role, QTL epistases were analyzed by four single-segment substitution lines with heading date QTLs in this paper. We found that QTLs of three positive effects and one negative effect generated 62.5% negative dual QTL epistatic effects and 57.7% positive triple QTL epistatic effects, forming the relationship "positive QTLs-negative one order interactions-positive two order interactions". In this way, the aggregation effect of QTLs was partially neutralized by the opposite epistatic effect sum. There also were two exceptions, QTL OsMADS50 and gene Hd3a-2 were always with consistent effect directions with their epistases, implying they could be employed in pyramiding breeding with different objectives. This study elucidated the mechanism of epistatic interactions among four QTLs and provided valuable genetic resources for improving heading date in rice.